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ECCENTRICALLY STIFFENED SHALLOW SHELLS 


OF DOUBLE CURVATURE* 

By John A. McElman 
Langley Research Center 

SUMMARY 

Equilibrium equations and boundary conditions are derived which govern the 
buckling and vibration of eccentrically stiffened shallow shells of double curvature. The 
equations are then solved for the case of simple- support boundary conditions and results 
are presented which illustrate the effects of eccentric stiffening on the dynamic and 
buckling characteristics of shells of positive and negative Gaussian curvature. Results 
show that eccentric stiffening can have a significant effect on the natural frequencies and 
buckling loads of such shells. 


INTRODUCTION 

Shells of double curvature have become common structural members in aerospace 
vehicles. The understanding of the effects of stiffening on the behavior of such shells has 
grown in importance with the increasing need for precision in the design of lightweight 
structures. The effects of stiffener eccentricities on the buckling and vibration character- 
istics of circular cylinders and flat plates have been studied analytically in references 1 
to 7. A typical stiffened cylinder is shown in figure 1. Experimental results of refer- 
ences 3 and 8 show that for some configurations a circular cylindrical shell stiffened with 
stringers attached only to its external surface can carry more than twice as much load in 
axial compression as its internally stiffened counterpart. 

In the present study, nonlinear equilibrium equations and boundary conditions are 
derived for eccentrically stiffened shells of double curvature. The types of doubly curved 
shells considered are shown in figure 2. These particular shell configurations are chosen 
because exact closed-form solutions can be obtained which exhibit typical eccentricity 
effects for shells of double curvature. The derivation is accomplished by defining non- 
linear strain-displacement relations for the shell and stiffeners. The potential energy of 


*The information presented herein is based in part upon a thesis offered in partial 
fulfillment of the requirements for the degree of Doctor of Philosophy in Engineering 
Mechanics, Virginia Polytechnic Institute, Blacksburg, Virginia, June 1966. 



the system is then formulated and the nonlinear equilibrium equations and boundary con- 
ditions are obtained by the application of the principle of minimum potential energy. The 
nonlinear equations are subsequently linearized to obtain equations which govern the 
small- amplitude vibrations of prestressed eccentrically stiffened shells. Closed-form 
solutions to these equations are presented for a membranelike prestress state and simple- 
support boundary conditions. If the natural frequency is assumed to be zero, the linear 
equations and the solution apply to buckling problems. 

Because of the large number of parameters involved, presentation of results of a 
general nature would be impractical. Results are included, however, for the free vibra- 
tion of specific shell configurations to demonstrate the effects of eccentric stiffening on 
vibration problems. Data are also presented for the buckling in axial compression and 
under hydrostatic pressure for specific shell configurations. These findings illustrate 
typical effects which exist in eccentrically stiffened shells of double curvature. 

SYMBOLS 

The units used for the physical quantities defined in this paper are given in both the 
U.S. Customary Units and the International System of Units (SI) (ref. 9). Factors relating 
these two systems are presented in appendix A. 

A cross-sectional area of stiffener 

a length of shell 

C^,C2,C^,C^ constants defined in appendix B 

Et^ 

D flexural stiffness of isotropic shell wall, — 

12(1 - M 2 ) 

d stringer spacing 

E Young's modulus 

e distance from shell middle surface to line on which N x acts 

e r ,e s nondimensional eccentricity parameters 

f frequency, 

2tt 

G shear modulus 
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I 

I moment of inertia of stiffener about an axis through its centroid parallel to 

the shell middle surface 

I Q moment of inertia of stiffener about middle surface of shell 

J torsional constant for stiffener 



u,v,w tangential displacements and normal displacement of shell middle surface 

in x-, y-, and z-directions, respectively 

u,v,w displacement amplitudes 

x,y,z rectangular Cartesian coordinates 

Z curvature parameter, ^-r(l - 


3 


distance from middle surface of shell to centroid of stiffener 


«,/3 


wavelength parameters 


r defined by equation (33) 

e x ,€y,yxy normal and shearing strains at shell middle surface 
e xT’ e yT’ y xyT total normal and shearing strains 

A,A 0 ,A r ,A s ,A rs eccentricity parameters defined by equations (34) 
Aj,A 2 ,A- 3,>.4 defined by equations (BIO) 


M 


n 


Poisson's ratio 


R 1 

curvature parameter, 
potential energy 


mass density 


O) 


circular frequency 


Subscripts: 


A prestress state 

B small changes away from prestress state 

i inertial 

L load 

r stiffening in y-direction (rings) 

s stiffening in x-direction (stringers) 
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sh 


shell 


T total 

A subscript preceded by a comma indicates partial differentiation with respect to 
the subscript. 


ANALYSIS 

In this section, the nonlinear equilibrium equations and boundary conditions for the 
shallow shells of double curvature illustrated in figure 2 are derived. Solutions are 
obtained for the case of simple- support boundary condition. 


Energy Formulation 

Strain-displacement relations .- The strain-displacement relations for the shells of 
double curvature are derived in reference 10 and are as follows: 

1 O w ^ 

e xT = u, x + 2 w ’x ± Rjj ~ zw ’xx 

e y T = v, y + f w, y + " zw ’yy | (i; 

y X yT = u ’y + v 'x - w >x w >y " 2zw ’xy 

J 

where u and v can be identified as the tangential displacements of the middle sur- 
face of the shallow shell and w can be regarded as the normal displacement. In equa- 
tions (1) and in subsequent equations, a double sign is used; the upper sign applies to 
shells with positive Gaussian curvature and the lower sign applies to shells with negative 
Gaussian curvature. 

Strain energy of isotropic shell .- The strain energy of the unstiffened thin-wall 
isotropic shell is 


nsh 2(l - p 2 p-t/2 *^0 ^0 


■ f‘ /2 f 2lRl f a f £ 2 T + 4 T + 2M, xT e yT + 


1 - M 
2 


rjy-rjdx dy dz (2) 


Substitution of equations (1) into equation (2) and integration with respect to z yields the 
following expression for shell strain energy: 
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( 3 ) 


n 


sh 



? + 2 ' 16 x £ y +i f Jir iy dxdy ) 

+ 2 ^ w ’xx w ’yy + ■ (‘'"'•ly)* 1 dy 


where 

e x = e xT| z=0 ’ e y = € yT| z=0 ’ ^xy = r xyT| z=0 (4) 

2+3 

and D = — -r is the flexural stiffness of the shell wall. 

12(l - M 2) 

Strain energy of stiffeners .- The strain energy of the stiffeners is derived by using 
the fact that the displacements in the shell and stiffeners are equal at the point of attach- 
ment. Integration with respect to the thickness coordinate yields terms whose signs 
depend upon whether the stiffeners are attached to the inner or the outer surface of the 
shell. The stiffeners are assumed to be in a state of uniaxial stress, and stiffener 
twisting is accounted for in an approximate manner. For configurations in which both 
rings and stringers are attached to the same surface of the shell, the effects of joints in 
the stiffener framework are ignored. 

The total strain energy of K stringers in the x-direction is written as 



where the first term within the parentheses is the strain energy of bending and extension 
in the stiffener, and the second term is the strain energy of twisting of the stiffener. This 
latter term results from assuming that the stiffener twists in a fashion such that its angle 
of twist is equal to the local angle of twist of the shell. The quantity dA s is an element 
of the cross-sectional area of the stiffener and the quantity G S J S is the twisting stiff- 
ness of the stringer section. After substitution from equations (1) and (4), the first term 
inside the parentheses can be written as 



Inspection of these terms reveals that the first integral inside the parentheses is the area 
A s of the stiffener cross section, the second integral is the first moment of the area 
z g Ag where z s is the distance from the middle surface of the isotropic shell (z = 0) to 
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the centroid of the stiffener cross section, and the third integral is the moment of inertia 
of the stringer I os about z = 0. Note that the centroid distance z s is positive for 
external stiffeners and negative for internal stiffeners. If the stiffener spacing d at the 
equator is sufficiently small, the effects of the stiffeners can be averaged or "smeared 
out," and the finite sum in equation (5) can be replaced by an integral. Equation (5) then 
becomes 




O \ VTgU g €\ 

s e x w >xx + w.y + - ^ w 


G S J 
- — — 

2 ’xyj 


dx dy 


( 6 ) 


The strain energy of stiffening in the y-direction can be derived in a similar manner and 
is written as follows: 


n r 



2z r A r 


€ y w >yy 


+ i or w, 


G r J r 


’yyj 


w 


’xyj 


dx dy 


( 7 ) 


where l is the ring spacing and the subscript r is used to denote ring properties. 

Potential energy of applied loads .- The only applied loads to be considered in this 
analysis are an external pressure p and an externally applied load resultant N x 
(positive in compression). The potential energy associated with these loads is 


p27rRq pa 

n L'J 0 J 0 pwdx 




(8) 


The quantity e is the distance from the middle surface of the isotropic shell (z = 0) to 
the line on which N x acts. 

Potential energy of in ertia loading .- If the stiffened shell is undergoing simple har- 
monic motion of circular frequency co (inplane inertia neglected), and if u, v, and w 
are amplitudes of such motion, the potential energy of inertia loading at maximum deflec- 
tion is 


n w = - i f 1 f Mw 2 w 2 dx dy (9) 

2 J Q Jq 

A s A r 

where M = p , t + p — — + p is the smeared-out mass per unit of surface area of the 

sn s d r j 

stiffened shell. The quantities p g ^, p g , and p r are the mass densities of the shell, 
x-direction stiffeners, and y-direction stiffeners, respectively. 

Nonlinear Equilibrium Equations and Boundary Conditions 

The total potential energy n T is the sum of the energies given by equations (3), 

(6), (7), (8), and (9): 
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( 10 ) 


= ^sh + Us + Hr + n L + n w 


Equation (10) may be written in terms of stress and moment resultants as 

i p 2irRj p& < . 

n T = 2 Jq Jq ^x e x + ^y e y + ^xy^xy " ^x w >xx " ^y w, yy + ^xy w ’xy " ■^yx w )xyj < ^ x dy 


+ I 0 2,,Rl IJV ' 5 M “ 2w2 ) dx dy + I 0 2 ’ Rl N *(“ - ew ’x) 
where the stress and moment resultants are defined as follows 
Et 


dy 


(ll) 


and 


N x = 


1 - /Ll 2 

E, 

+ 


E S A S 


Ny = 


d 

Et 


1 - ii* 

E r A r 
+ l 


u ’x * 5J + \ W -I + <* ( v ’y + iq- + \ w ’?) 
( U ’x ± ^ t | W ’x- 5sW ’xx) 

v ’y * % * I w ’y * t *( u ’ x * r| + 1 

\ 


v ’y + ^ + l w ’?- 5 >- w ’yy. 


Nxy Gt^l, 


,y + v, x + w, x w,y 


) 


My = 


My = - 


Vs , ri ^s E s^s 
d ,xx ■ d 


E r!r 


M 


xy 


D ( w ’xx + MW,yyj + 

D (w, yy+ 

; ) w, xy 


w 1 9 

U ’X ± Ri + 2 W ’X Z S W ’XX 


> (12) 


z r E r A r 


wl9- 

f^i v, y + rT + 2 w ’y ■ z r w ’yy 


M. 


yx 


■(¥■ 


G S J S 


G r J r 


w. 


i / ,xy 


Is _ Iqs " ^s^s 


*r *or “ ^r^r 
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The stress resultants N x and Ny are positive in tension. The nonlinear equilibrium 
equations and boundary conditions are obtained from equation (11) by application of the 
principle of minimum potential energy (611 = 0) and the fundamental lemma of the calculus 
of variations. The equations so obtained are 

N x,x + N xy,y = 0 


N y,y + N xy,x “ 0 


-NU 


+ M - M 

xy,xy yx,xy 


Nx Ny 

M y,yy * r£ + rJ 


- N X W >XX - NyW,yy - 2N x yW, X y - MO)^ + P = 0 


(13) 


A set of boundary conditions to be satisfied at each end of the shell (x = 0,a) are 

Mv -V — (Mwt 


Myx,yj + N x w, x + NjjyW,y - 0 

(14a) 

or w = 0 

(14b) 

M x + N x e = 0 

(15a) 

or w, x = 0 

(15b) 

N x + N x = 0 

(16a) 

or u = 0 

(16b) 

Nxy = 0 

(17a) 

or v = 0 

(17b) 


The natural boundary conditions are given by the expressions in equations (14a), (15a), 

(16a), and (17a) and the geometric boundary conditions are given by equations (14b), (15b), 
(16b), and (17b). The condition in equation (14a) requires that a shear resultant comparable 
to the Kirchhoff shear vanish and hence is a free-edge boundary condition. The three 
natural boundary conditions in equations (15a), (16a), and (17a) correspond to conditions 
in which the edge moment resultant, the extensional stress resultant, and the shear stress 
resultant, respectively, vanish. 


Homogeneous Equations Governing a Prestressed Vibrating Shell 

In this section, the nonlinear equilibrium equations (eqs. (13)) are used to obtain 
linear equations which govern the small- amplitude vibration of a prestressed eccentrically 
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stiffened shallow shell of double curvature. The deformations u, v, and w associated 
with the vibration of a prestressed shell are divided into two parts as follows: 

u = u A + u b> v = V A + V B’ w = W A + W B (! 8 ) 


The first part, denoted by the subscript A, is assumed to be an axisymmetric static pre- 
stress deformation which occurs prior to the excitation of one of the natural frequencies. 
The second part, denoted by the subscript B, is a small additional deformation which 
occurs as a result of the excitation. Since equations (13) are equilibrium equations for 
the system, the displacements denoted by subscript A as well as the sum of the two dis- 
placements denoted by subscripts A and B must satisfy these equations. After sub- 
stitution of the axisymmetric (subscript A) displacements, equations (13) become 


N xA,x = 0 
N xyA,x = 0 

N xA NyA 

- M xA,xx ± iff + if - N xA w A,xx + P = 0 




(19) 


N. 


xA 


N, 


yA 


now 


Et 

11 

l-M 2 

A, 

Et 

™A 

lV 

Rl 

= Gtv A,x 


W A . 1 ...2 


/w. 




w. 


Ri, 


E S A S , 


w A i 2 - 1 

l U A,x ± + 2 W A,x ' z s w A,xx) 


A 12 ] 

iI T TA,x^ + 2 W A,xi 


E r A r r Wp\ 

i IrJ 


> ( 20 ) 


M 


xA 


E S I S 

Dw A,xx + “ d W A,xx 


ZsE s A s 


W A 1 2 - \ 

, u A,x * Rj + 2 w A,x - z s w A,xxJ 


J 


A set of appropriate boundary conditions is found from equations (14) to (17) to be 


M xA,x + N xA w A,x = 0 or W A = 0 
M xA + *x e = 0 or w A ,x = 0 
N xA + N x = 0 or u A = 0 

N xyA = 0 ° r V A =0 


( 21 ) 
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A solution to equations (19) satisfying the conditions (21) completely describes the pre- 
stressed state. 


The equilibrium equations governing the additional deformations (subscript B) are 
obtained by substituting equations (18) into equations (13). If only linear terms in the 
additional deformations are retained and equations (19) are considered, equations (13) 
become 


N xB,x + N xyB,y " 0 
N yB,y + N xyB,x = 0 


-M. 


xB,xx + M xyB,xy " M yxB,xy - M yB,yy ± R 2 + Ri 


- M 


M. 


N xB . N yB 


N aW„ _ - N a w t 


N—W ; 


xA B,xx yA B,yy ^xB w A,xx 
and the boundary conditions become 


- MoiTfg = 0 


( 22 ) 


M. 


xB,x ' ( M xyB,y " M yxB,y) + N xA w B,x + N xB w A,x = 0 or W B = 


M xB = 0 or w B,x = 0 


N xB = 0 or U B = 0 


N xyB ~ 0 or V B = 0 


(23) 


where 


N. 


Et 


xB _ 1 - M 2 


W I 

U B,x ± ^ T w A,x w B,x 


w T 


+ w a v w n T + M v„ + 


B 

B ’ y " Rly 


E s A s f W B - ' 

+ “ ~ \ U B,x ± R^ + W A,x w B,x " z s w B,xx, 


(24a) 


N. 


Et 


yB 


i - ^ 


w T 


' ra / Wg i 

^B,y + ^ + H u B,x * ^ + W A,x w B,x 


E r A r 


( VB - 


y + Ri " Z r W B,yy) ( 24b ) 


N xyB = Gt ( u B,y + v B,x + w A,x w B,y) 


(24c) 
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(24d) 


M xB “ ‘ 


|°(»B,x* + -B, yy ) ^ ^ »B,ss * ^ ( u B,x * i 


+ w A,x w B,x - z s w B,xx 


M. 


yB 


D ( w r „„ + MW 


E r I r 
+ — =-^- w 


z r E r A r 


w 


V„ „ + — ^ - Z^W 


B,yy ^ B,xxJ z B,yy t \’B,y T R x ~r" B ,yy 


M. 


xyB 


= (ot! + 2s£s) 

\ 6 d / B,xy 


(24 e) 


(24f) 


“yxB'-^^Kxy 


(24g) 


The homogeneous equations (22) and the homogeneous boundary conditions (23) represent 
an eigenvalue problem which governs the natural frequencies of a prestressed eccentri- 
cally stiffened shell of double curvature. The coefficients in the equations are determined 
by considering solutions to the axisymmetric prestress problem described by equa- 
tions (19) and the associated boundary conditions (eqs. (21). 


Solution for Prestressed Vibrating Shell 


Equations (22) have variable coefficients and would be quite difficult to solve in 
most instances. If, however, an assumption analogous to that made in classical buckling 
theory is made (i.e., that the lateral prestress deformation w^ is constant prior to the 
excitation), the solution is greatly simplified. The implications of this assumption are 
given in reference 6, and the exact solution to equations (19) and (20) is presented in 
appendix B. 


For w^ constant and no applied shear, the prestress equations are 


N x ^ = Constant I 


= 0 


NxyA = Constant 


> 


(25) 



From the boundary conditions (eqs. 21)) 


N xA " " N x 


Substituting equation (26) into the third of equations (25) yields 

N y A = "(p R 1 * *V) 


R 


1 


where v = — - is a curvature parameter. 

R2 

Also for w^ constant, equations (22) become 


N xB,x + N xyB,y " 0 
N yB,y + N xyB,x = 0 


N xB n v b 

_M xB,xx + M xyB,xy “ M yxB,xy ' M yB,yy ± ^ + 


~ N xA w B,xx - N yA w B,yy ’ Mw W B = 0 




where now 


N. 


Et 


xB " 1 - M 2 


Wg / w B 

Ut, _ ± — — + U Vtj „ + 

B,X R 2 H B,y Rl 


E s A s l W B - 

+ — r— Ut, „ ± - z e w, 


d rB,x-R^ s B, 


XX/ 


N. 


Et 


y B " 1 - 


W B / W B ' 

J,y + rT +m ( Ub ,x ± ^ / 


E r A r / w B 
+ — ( v B,y + Rj- z rWB,yyl 


N xyB - Gt ( u B,y + v B,x) 


M. 


xB 


. EoIq ZcEoAo / w R 

D ( W B,xx + ^B,yy) + — w B rXX ‘ ~ 5“ K.x * rJ * z s»B, 


.XX/ 


M 


yB 


E r Ir z r E r A r / W B 

D ( W B,yy + < lw B, X x) + — "B,yy - — — ( v B,y + R^ ‘ 


B V'B.yy 


(26) 

(27) 

(28) 

(29a) 

(29b) 

(29c) 

(29d) 

(29e) 
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,, _ (Gt 3 G s J s\ 

M xyB ( 6 + d / B,xy 


M. 


yxB 


(Gt 3 G r J r\ 
f» + l / W B,xy 


(29f) 

(29g) 


If the origin of the coordinate system is taken at one edge of the shell, the simple- 
support boundary conditions to be satisfied are 


? B(0,y) xB(0,y) ~ V B(0,y) " W xB(0,y) 

7 B(a,y) _ M xB(a,y) = V B(a,y) = N xB(a,y) 


(30) 


Expressions for the displacements Ug, Vg, and Wg which satisfy these boundary con- 
ditions are given as 


Ug = u cos cos ny 




Hi 


v R = v sin sin — 
b a Rj 


w = w sin ES cos 2 - 
B a Ri 

v 


(31) 


where m is the number of axial half waves and n is the number of circumferential 
full waves. After substitution of equations (31) into equations (28) and nondimensionaliza- 
tion, the following equation is obtained: 


[An] 

- ju 2 ) + IX ± v ± Sv(l - jli 2 )J 
[e r R(l - m 2 ) n 2 ± fiv + 1 + (l - ju 2 )r] 


e s S c?(l - (i 2 ) + ju ± v ± Si>(l - ju. 2 )] Je r R(l - ju 2 )n 2 ± [iv + 1 + (l - ju 2 )r] 


-£ + 5(i -m 2 ) + ^/ 


■M 


1 + fl 


1 + r(i - n 2 ) + i^Ji 

_ 2 


L 2 / 3 2 J 



— 


~ - 


w 


0 


u 

= 

0 


V 


0 

J 




__ 


( 32 ) 
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where 


_ Doft(l - /i 2 )(l + ft 2 )"* ^yA n ^ (-*• ~ M^) ^xA 0 ^ " ^ 2 ) 

EtR 2 Et Et 

/G S J S | G r J r y 2n 2( 1 _ u 2) E r I r (l - M 2 )n 4 E s I s (l - m 2 )^ . 2 w 2) _ 

\ d l ! EtR 2 l EtR 2 EtdR 2 ^ M 

- e 2 n 4 (l - h 2 )r - 2e r n 2 (l - m 2 )r t 2a 2 z,e s (l - m 2 )s - (l ± 2[iu + u 2 ) 


-( o\ -/ o\ o Ma) 2 R?(l - M 2 ) 

- R(l - M 2 ) - s(l - \i 2 )u 2 + 1 

Et 


and the following nondimensional parameters have been defined: 


e 




m7jRi 

“ = 0 = 


na 

m7rRi’ 


E r A 
EtZ ’ 


S = 


E s A s 

Etd 


For a nontrivial solution to exist, the determinant of the coefficients of u, v, and w 
must equal zero. After some manipulation, the nondimensional characteristic equation 
becomes 


Mo> 2 a 4 

7T 4 D 




r = m 4 


(l + f> 2 )‘ 


+ m 


Dd 


E sh + ^ ^s d s + ^r d r^ 2 + ®r*r 


12Z 2 f A o + SA s + RA r + RSA 

7T 4 \ A 


cLD 


rs 


ZD 


DZ 


(33) 


where 


A 0 = (l ± up 2 ) 2 


(34a) 

A s = 1 ± i//3 2 (2 ± up 2 ) 

1 + 2e s a 2 (y3 2 - jn)(l ± up 2 ) + e 2 a 4 (l + p 2 f 

(34b) 

A r = 1 ± vfi 2 (2 ± up 2 ) 

1 + 2e r n 2 (l - np 2 )(l ± up 2 ) + e 2 n 4 (l + p 2 ) 2 

(34c) 


15 



A rs = egn 2 ^]^ 1 - n 2 )p 2 + 2(1 + m)] + e 2 n 4 £(l - n 2 ) + 2/3 2 (l + m)] + 2e r e g n 4 (l + /i) 2 

+ 2(l - /i 2 )n 2 (l ± t>/3 2 )e r + 2(1 - p. 2 )n 2 (l ± i/)3 2 )e s + (l - /j. 2 ) jl ± t»/3 2 (2 ± t»/3 2 )] (34d) 

A = (j3 2 + l) 2 + (R + S)(l + m) 2 0 2 + (l - m 2 )[s + R/3 4 + 2/3 2 RS(l + /i)] (34e) 


and the nondimensional parameter 


z 2 - a 4 (l - u 2 ) 
R 2 t 2 


has been defined. 


Equation (33) is a closed-form expression which gives the natural frequencies for a 
prestressed eccentrically stiffened shell of double curvature. If the natural frequency 
co is set equal to zero in equation (33), a stability equation results which may be mini- 
mized to obtain buckling loads or buckling coefficients for a variety of loadings. 

In equation (33), the effect of eccentric stiffening is reflected by the terms con- 
taining the quantities e r and e g . These quantities are positive for stiffening attached 
to the external surface of the shell and negative for stiffening attached to the internal sur- 
face. The term (/3 2 - jli) in equation (34b) and the term (l - pt/3 2 ) in equation (34c) may 
also change sign, depending upon the shell geometry and its deflected shape. It should 
also be noted that every term that is first degree in e r or e s is modified by the term 
(l ± vfi 2 ), an indication that the type and magnitude of the shell curvature can greatly 
influence the effects exhibited by eccentric stiffening. 


RESULTS AND DISCUSSION 


The results presented include the vibration of shells in the absence of prestress, 

the buckling of shells in axial compression, and the buckling of shells under a hydro- 

_ PRi 

static pressure loading with N x = — . The prestress load resultants for axial compres- 
sion are: 


N xA = - N x 


N . = ±N v\ 
yA x 1 

and the prestress load resultants for hydrostatic pressure are: 


N. 


xA 




(35) 


( 36 ) 
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The following characteristic equations for each of these cases are found from equa- 
tion (33) by using equations (35) and (36): 


M 2 4 

Vibration (no prestress): — j T 

7/ D 


Buckling (axial compression): 


N x a 2 


Dir 2 m(l T 1^/3 2 ) 


Buckling (hydrostatic pressure): 


pRja^ 

Dir 2 


( 37 ) 


(38) 


(39) 


Numerical results obtained from these equations are discussed in the following sections. 


Vibration 

All the vibration results presented for shells of double curvature were obtained 
from equation (37) for a shell whose stiffener configuration is shown in figure 3. The 
natural frequencies in the absence of prestress for a stringer- stiffened shell of positive 
Gaussian curvature are shown in figure 4. The natural frequencies for external stiff- 
ening are in general higher than those for internal stiffening, the difference being about 
35 percent for m = 2 and n = 8. A crossover does occur at low values of n due to 
the fact that the term (/3 2 - [i) has changed sign in equation (34b). The corresponding 
frequencies for a cylinder {v = 0) are shown in figure 1 of reference 7. The shell of 
positive Gaussian curvature is stiffer and as a result, it exhibits higher frequencies than 
the corresponding cylinder. 

Figure 5 shows the natural frequencies in the absence of prestress for a stringer 
stiffened shell of negative Gaussian curvature (u = -0.25). The curves for external and 
internal stiffeners cross twice for each m. The first crossing occurs when the term 
(/3 2 - ju.) changes sign, whereas the second crossing occurs when the term (l - v/3 2 ) 
changes sign. This trend makes it difficult to determine which type of stiffening yields 
the highest frequency for a given m and n. The figure does illustrate, however, that 
stiffening eccentricities can significantly alter the natural frequencies of a stiffened shell 
of this type. External stiffening yields a natural frequency that is approximately 
50 percent higher than that for internal stiffening for m = 2 and n = 4, and internal 
stiffening gives approximately a 50-percent higher frequency at m = 2 and n = 8. The 
analysis of the shell of negative Gaussian curvature is thus more complicated than that of 
either the corresponding shell of positive Gaussian curvature or the cylinder (ref. 7). 

It is necessary to note that the results depend on the shell geometry since /3 is a 
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function of the shell length and radius as well as mode shape and they also depend on the 
ratio of principal curvatures. The crossings in the figure can always be predicted, how- 
ever, by setting the terms (/3 2 - p) and (l - v($) equal to zero. 

Figures 6 and 7 show the natural frequencies with no prestress of ring stiffened 
shells with v = 0.25 and v = -0.25, respectively. As in the case of ring stiffened 
cylinders, the shells of double curvature considered herein exhibit relatively small eccen- 
tricity effects when stiffened by rings alone, the maximum eccentricity effect being about 
10 percent. The points where the curves cross in the figure can be predicted by consid- 
ering the terms (l - p/3 2 ) for figure 6 and (l - p0 2 ) and (l - i//3 2 ) for figure 7. 

Buckling 

The buckling results obtained from the present analysis are for a shell with a stiff- 
ener configuration like that shown in figure 3 loaded either by end loads or hydrostatic 
pressure. Tables I and II give the pressures and the compressive end loads for buckling 
of shells of positive ( v = 0.25) and negative {v = -0.25) Gaussian curvatures, respectively. 
The first column of results in the tables gives the critical axial load resultants for a shell 
loaded in axial compression, and the second column of results gives the critical values of 
pressure for a shell loaded by hydrostatic pressure. The effects of eccentric stiffening 
are evident from a consideration of the tables; a few important facts should be mentioned, 
however. The eccentricity effects are not as large for the shell loaded in axial compres- 
sion as they are for the corresponding cylinder (ref. 6). It is also evident that stringers 
are more effective stiffeners than rings under axial compressive loadings for the shell of 
positive Gaussian curvature but rings are more effective stiffeners than stringers under 
this type of loading for a shell of negative Gaussian curvature. Rings are more effective 
than stringers under hydrostatic pressure loading for both types of shells, as is true for 
cylinders (ref. 6). 

The data of table I show that external stringers are more effective (approximately 
20 percent) than internal stringers under axial compressive loading for a shell of positive 
Gaussian curvature with no ring stiffeners. However, internal stringers are more effec- 
tive under this loading for a shell of negative Gaussian curvature, as shown by the data 
of table n. 

Examination of equations (35) and (33) shows that a shell of positive Gaussian curva- 
ture can buckle with an applied tensile edge load because of the buildup of compressive 
hoop stresses. Table III gives the tensile buckling loads for such a shell ( v = 0.25). It 
is evident from the table that rings provide the most effective type of stiffening for this 
case with internal rings giving a 28-percent higher buckling load than external rings. 
Combining both types of stiffening does not greatly increase the buckling load and, in 
fact, table III shows that external rings and internal stringers give a lower buckling load 
than internal rings alone. 
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CONCLUDING REMARKS 


An analysis is made of the buckling and vibration of eccentrically stiffened shallow 
shells of double curvature. An expression is presented relating the natural frequencies to 
prestress terms and to a variety of nondimensional shell and stiffening parameters. 

This expression may be used to determine the vibration and buckling characteristics of 
particular shell- stiffener configurations or to perform parametric studies to optimize 
the stiffening configuration for a particular application. 

All the results presented herein illustrate the complicated behavior of eccentri- 
cally stiffened shells. It is impossible to make generalizations regarding the dynamic or 
buckling behavior of such shells. The eccentricity effects depend on the type of loading, 
the configuration, and the physical properties of both the shell and the stiffening. As a 
result, each particular shell- stiffening configuration must be thoroughly analyzed to 
determine these effects. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 30, 1966, 

124-11-06-04-23. 


APPENDIX A 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 


The International System of Units (SI) was adopted by the Eleventh General 
Conference on Weights and Measures held in Paris, October 1960 in Resolution No. 12 
(ref. 9). Conversion factors for the units used herein are given in the following table: 


Physical quantity 

U.S. Customary 
Unit 

Conversion 

factor 

(*) 

SI Unit 

Force 

lbf 

4.448 

newtons (N) 

Frequency 

cps 

1 

hertz (Hz) 

Length 

in. 

0.0254 

meters (m) 

Stress and pressure 

lbf/in 2 

6.895 X 10 3 

newtons/ meter2 (N/m2) 

Unit loading 

lbf/in. 

175.1 

newtons/meter (N/m) 


^Multiply value given in U.S. Customary Unit by conversion factor to obtain 
equivalent value in SI Unit. 

Prefixes to indicate multiple of units are as follows: 


Prefix 

Multiple 

centi (c) 

1 

O 

iH 

mega (M) 

10 6 

giga (G) 

10 9 
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APPENDIX B 


PRESTRESS DEFORMATIONS 


The solutions presented in this report are based on the assumption of constant lat- 
eral displacement w A prior to buckling. The implications of this assumption for a 
cylindrical shell are discussed in reference 6. The discussion in reference 6 applies 
directly to the shells of double curvature except that the added curvature can possibly 
increase the importance of considering exact prestress deformations rather than the 
simple assumption of constant lateral displacement. Consideration of exact prestress 
deformations would have necessitated an approximate solution to equations (22) for 
buckling and vibrating. An exact prestress solution can be obtained, however, and is 
presented as follows. 


The equations which govern axisymmetric prestress deformations are 

n xA,x * 0 


N xyA,x 0 


(Bl) 

(B2) 


where 


M + 

-M xA,xx ± R 2 + R! 


N xA w A,xx + P = 0 


(B3) 


E„A e 


xA 1 _ m 2( € xA + Me yA) + a ( e xA “ z s w A,xx) 

v=7^( £ y A + ^) + ^M 


M. 


xA 


Dw 


^s^s 


A, xx ’ r d W A,xx 


z s E sAs 




e xA - z s w A ,xx) 


(B4) 

(B5) 

(B6) 


Equation (Bl) implies that N xA = Constant. This constant is denoted by -N x , where 
N x is an applied compressive end load. In addition, equation (B2) implies that 
N A = Constant = 0 with no applied shear load. For axisymmetric deformations, 

J W A 

e yA = thus equation (B4) can be solved for as follows: 


E s-^s^s 


w Etfx W A jT 

d A, xx ! _ u 2 Ri x 

Et E S A S — 


: xA 


(B7) 


1 - M 2 
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Substitution of equation (B7) into equations (B4) to (B6), and subsequent substitution 
of equations (B4) to (B6) into equation (B3) yields 


where 


w A 

A,xxxx 


c 2 C 3 

+ Cj W A,xx + cY W A 



Cl - P + 


E sls 


[l + s(l- m 2 )] 


^s E s A s 


C 2 - 22 R^ tM + N x[j + S(l - M 2 )] 


_ Et n 

x, 9 . li 


C 3 - 2 

Rl 2 


1 + R + S + RS 


(l - M 2 ) 


C4 = ( p T ^IC 1 + s ( 1 - <* 2 )] 


Rl 


A solution to equation (B8) is of the following form: 


(B8) 


w. 


XjX 


AoX AoX XaX 

' ■ X 3 e 6 + K 4 e 4 


_4 

C 3 


(B9) 


where the K's are constants to be determined from the boundary conditions and the 
X's are defined as follows: 



(BlOa) 



(BlOb) 


x 3 = 



(B 10c) 
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(BlOd) 


Equation (B9) can be used in conjunction with equations (22) to determine the effect 
of prestress deformations on the buckling and vibrating characteristics of doubly curved 
shallow shells. 
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TABLE I.- BUCKLING RESULTS FOR SHELL OF POSITIVE GAUSSIAN CURVATURE 


Shell properties 

v = 0.25 

l = d = 1.0 in. (2.54 cm) 

E = E s = E r = 10.5 x 10 6 psi (72.4 GN/m 2 ) 
a = 23.75 in. (60.3 cm) 

R 1 = 9 .55 in. (24.3 cm) 

/z = 0.3 


Configuration 

Axial compressive 
load, N x 

Hydrostatic pressure 
load, p 

lbf/in. 

MN/m 

psi 

MN/m 2 

No stringers 
Rings external 

744 

0.130 

156 

1.076 

No stringers 
Rings internal 

717 

0.126 

150 

1.034 

No rings 

Stringers external 

7,913 

1.385 

24.5 

0.169 

No rings 

Stringers internal 

6,664 

1.166 

22.1 

0.152 

Rings external 
Stringers internal 

10,144 

1.776 

374 

2.580 

Rings external 
Stringers external 

11,778 

2.063 

413 

2.849 

Rings internal 
Stringers internal 

7,928 

1.389 

415 

2.860 

Rings internal 
Stringers external 

9,183 

1.607 

462 

3.187 



TABLE II.- BUCKLING RESULTS FOR SHELL OF NEGATIVE GAUSSIAN CURVATURE 


Shell properties 

v = -0.25 

l = d = 1.0 in. (2.54 cm) 

E = E s = E r = 10.5 x 10 6 cm (72.4 GN/m 2 ) 
a = 23.75 in. (60.3 cm) 

R-j^ = 9.55 in. (24.3 cm) 
ji = 0.3 


Configuration 

Axial compressive 
load, N x 

Hydrostatic pressure 
load, p 

lbf/in. 

MN/m 

psi 

MN/m 2 

No stringers 
Rings external 

744 

0.130 

69 

0.476 

No stringers 
Rings internal 

697 

0.122 

63 

0.434 

No rings 

Stringers external 

282 

0.049 

14.5 

0.100 

No rings 

Stringers internal 

340 

0.060 

18.5 

0.128 

Rings external 
Stringers internal 

2,040 

0.354 

78 

0.538 

Rings external 
Stringers external 

1,913 

0.335 

70 

0.483 

Rings internal 
Stringers internal 

1,770 

0.310 

84 

0.579 

Rings internal 
Stringers external 

1,842 

0.322 

68 

0.469 
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TABLE III.- BUCKLING RESULTS FOR SHELL OF POSITIVE GAUSSIAN CURVATURE 

LOADED IN TENSION 

Shell properties 

v - 0.25 

I = d = 1.0 in. (2.54 cm) 

E = E s = E r = 10.5 X 10 6 cm (72.4 GN/m 2 ) 
a = 23.75 in. (60.3 cm) 

Rj = 9.55 in. (24.3 cm) 

II = 0.3 


Configuration 

Axial tensile load, 

N x 


lbf/in. 

MN/m 

No stringers 
Rings external 

13,441 

2.354 

No stringers 
Rings internal 

17,200 

3.010 

No rings 

Stringers external 

833 

0.146 

No rings 

Stringers internal 

751 

0.132 

Rings external 
Stringers internal 

16,617 

2.910 

Rings external 
Stringers external 

17,788 

3.115 

Rings internal 
Stringers internal 

19,408 

3.400 

Rings internal 
Stringers external 

20,396 

3.570 
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Rings 


Figure 3.- Stiffening configuration. All dimensions are in inches (millimeters). 


30 



PROPERTIES 


d = 1.0 in. (2.54 on) 

E = E g - 10*5 x 10^ psi (72.4 GN/m 2 ) 

a = 23.75 in. ( 60.3 cm) 

R 1 = 9.55 in. (24.3 cm) 

n = 0.3 

Stiffener configuration shown in figure 3 



— — External 
— — Internal 






f, cps. 

(Hz) 



n 


Figure 6.- Natural frequencies of ring- stiffened shell of positive Gaussian curvature (v 



12 3 4 5 


n 

Figure 7.- Natural frequencies of ring- stiffened shell of negative Gaussian curvature (v = -0.23). 
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